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ABSTRACT 

We probe the relationship between star formation rate (SFR) and radio synchrotron luminosity 
in galaxies at < z < 2 within the northern Spitzer Wide-area Infrared Extragalactic survey 
(SWIRE) fields, in order to investigate some of the assumptions that go into calculating the 
star formation history of the Universe from deep radio observations. We present new 610- 
MHz Giant Metrewave Radio Telescope (GMRT) observations of the European Large- Area 
ISO Survey (ELAIS)-N2 field, and using this data, along with previous GMRT surveys carried 
out in the ELAIS-Nl and Lockman Hole regions, we construct a sample of galaxies which 
have redshift and SFR information available from the SWIRE survey. We test whether the 
local relationship between SFR and radio luminosity is applicable to z = 2 galaxies, and look 
for evolution in this relationship with both redshift and SFR in order to examine whether the 
physical processes which lead to synchrotron radiation have remained the same since the peak 
of star formation in the Universe. We find that the local calibration between radio luminos- 
ity and star formation can be successfully applied to radio- selected high-redshift, high-SFR 
galaxies, although we identify a small number of sources where this may not be the case; these 
sources show evidence for inaccurate estimations of their SFR, but there may also be some 
contribution from physical effects such as the recent onset of starburst activity, or suppression 
of the radio luminosity within these galaxies. 
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1 INTRODUCTION 

In recent years there has been much interest in constrain- 
ing t h e star formation hi s tory of the Un iverse (e.g. iLillv et al 



199^; 'Madau et al." '1996'; "Mobasher et al. 1999; Haars ma et al 
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Magnelli et al. 2009). Estimates of the variation in the co-moving 
star formation rate density (SFRD) of the Universe with redshift 
can be made from a number of multi-wavelength tracers, with 
the most common on es being in the optical and infrared regimes 
(see lKennicut3ll998l . for a review of star formation indicators). 
The se estimates can differ by up to an order of magnitude (e.g. 
iHop kins et al. 2003), which is in part due to the significant uncer- 
tainty in the correction factor that needs to be applied for the effects 
of interstellar dust grains. These grains absorb short- wavelength 
ultra-violet (UV) and optical radiation, heat up to temperatures of a 
few tens to hundreds of kelvin, and re-emit thermally with a mod- 
ified black-body spectrum that peaks in the infrared, thus having a 
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significant effect on the optical and infrared spectral energy distri- 
bution (SED) of a galaxy. 

There have been a few attempts to estimate the SFRD of 
the Universe from radio observations (e.g. 'Mobasher et al.' n^99l : 
[Haarsma et al. 2000; Seymour etal. 2008; Smolcic et al. 20o3). 
These studies have a significant advantage over optical and in- 
frared estimates, in that radio observations are unaffected by the 
presence of dust grains, and do not require the corrections that are 
needed at other wavelengths. However, there are a number of dis- 
advantages to radio studies, the principal one being that very deep 
observations are required in order to reach the regime where star- 
for ming galaxies begin t o dominate (around 100 /iJy at 1.4 GHz, 
e.g. lSevmour et al .l2008h . This disadvantage will be overcome with 
the next generation of radio telescopes such as the Low Frequency 
Array (LOFAR) and the Square Kilometre Array (SKA), which are 
expected to routinely reach observation depths where the radio sky 
is dominated by star-forming galaxies. 

The non-thermal radio emission from star-forming galaxies is 
much greater than the thermal emission at frequencies below a few 
GHz, and is predominantly made up of synchrotron radiation, with 
supernova remnants putting energy into a population of electrons 
which are accelerated through the galactic magnetic field (see e.g. 
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ICondonlfTooi) . The details of this process are still a subject for 
debate, with two principal alternatives; that g alaxies are o ptically 
thick to cosmic rays and act as 'ca lorimeters' ( Volk 1989), or that 
galax ies are optically thin (IChi & W olfendale 1990; Helou & Bicavl 
Il993h . and the cosmic rays can diffuse out of the host galaxy. In the 
former case, the radio flux from a g alaxy will be esse ntially inde- 
pendent of the magnetic field B fe.g. lXhompson et aD [2006\ since 
all of the energy from the cosmic rays will be radiated within the 
galaxy. In the latter case only some fraction of the total energy will 
be radiated within the galaxy, and variations in the 5-field strength 
of galaxies over time may show up as a variation in the confine- 
ment of cosmic rays, and an evolution in the relationship between 
star formation rate (SFR) and synchrotron luminosity with redshift. 
Star-forming galaxies with magnetic field strengths of a few 
have been indirectl y detected through Faraday rotation mea sure- 
ments out to z 2 dlCronberg et allllOOSl : iBernet et al.ll2008h . but 
the rate and processes which lead to evolution in B-field strength 
with time are currently unknown. 

Attempts to estimate the SFRD of the Universe from radio ob- 
servations currently rely upon the untested assumpti on that local re- 
lationships between radio luminosity and SFR (e.g. Condon & Yin 
[l99Q : Bell 2003) can be applied successfully at higher redshift, and 
to sources undergoing much more extreme starbursts than are seen 
in the local Universe. In this work we test these assumptions, using 
new radio surveys of the three northern Spitzer W ide-area Infrared 
Extragalactic survey fSWIRE: Lonsdal e et al.1200 3) fields. SWIRE 
is a 49 deg^ study of six regions of sky that have particularly 
low values of background cirru s emission, using the Spitzer Space 
Telescope (I Werner et al .1120041) . Observations have been made in 
seven Spitzer wavebands with th e Multiband Imaging Array for 
Spitzer (MIP S: iRieke et all 12004) and the Infrared Array Camera 
(IRAC; Fazio et al.ll2004l) between 3.6 ^m and 160 jim. The three 
SWIRE fields studied in this work are in the northern hemisphere 
- the European Large Area ISO Survey-North 1 and -North 2 re- 
gions (ELAIS-Nl and ELAIS-N2), and the Lockman Hole. A great 
deal of complementary data has been taken on the SWIRE fields; 
a band-merged catalogue con taining over one n iillion sources has 
been created (Data Release 2: ISurace et al.ll2005|). and a photomet- 
ric r edshift catalogue has been generated feow an-Robinson et al.l 
I2OO8.) from the multi-wavelength data. All three northern regions 
are cov ered by the Faint Im ages of the Radio Sky at Twenty-cm 
(FIRST: iBecker et al.ll 19951) and NRAO VLA Sky Survey (NVSS; 
ICondon et al.lll998l) 1.4-GHz surveys; however the relatively shal- 
low sensitivity levels of these observations (0.75 and 2.25 mJy re- 
spectively) means that that these surveys are dominated by classi- 
cal radio galaxies, and contain few star-forming galaxies. We have 
been carrying out a series of deeper radio surveys of the northern 
SWIRE fields, in order to extend the multi-wavelength coverage 
currently available on star-forming galaxies into radio frequencies. 

In Section|2]we present new observations of the ELAIS-N2 re- 
gion, tak en at 610 MHz with the Giant Metrewave Radio Telescope 
(GMRT: AnanthakrishnanI [2005b . We combine these observations 
with our previous surveys of the ELAIS-Nl and Lockman Hole re- 
gions jGarn et al.ll2008al lbl) in Section[3] and create a sample of 510 
star-forming galaxies with photometric redshifts that are detected 
in the infrared and radio. In Section |4] we look at the infrared / ra- 
dio flux density ratio q'^^, and use deviations away from the typical 
values seen for star-forming galaxies to identify sources that may 
contain a significant excess of radio emission due to the presence 
of an active galactic nucleus (AGN). We compare our sample t o the 
local r elationship between radio luminosity and SFR given bv lBelll 
(l2003h . and find that the majority of sources are well described by 



this relationship, but that there are a small fraction of sources where 
the SFR would be significantly under-estimated by using the radio 
luminosity. In Section [5] we test the assumption that the local rela- 
tionship between SFR and radio luminosity can be used to calculate 
the SFRD of the Universe from galaxies at high redshift (z = 2), and 
from galaxies undergoing much more vigorous starbursts than are 
seen in the local Universe. 

Throughout this work a flat cosmology with the best-fitting 
parameters from the Wilkinson Microwave Anisotropy Probe 
(WMAP) fi ve-year data o^p.^ = 0.74 and Hq = 72 km s-^ Mpc^^ 
is assumed (iDunklev et al.ll2009h . 



2 RADIO OBSERVATIONS 
2.1 The ELAIS-N2 region 

The ELAIS-N2 region was observed in 2006 July 15 and 17 with 
the GMRT. Observations were made in two 16-MHz sidebands cen- 
tred on 610 MHz, each split into 128 spectral channels, with a 
16.9 s integration time. Thirteen pointings were observed, centred 
on 16^36"^48^ +41°01'45'' (J2000 epoch) and spaced by 36 ar- 
cmin in a hexagonal pattern, in a series of short interleaved scans in 
order to maximise the uv coverage. The typical time spent on each 
pointing was 22 min. A nearby phase calibrator, J1613+342, was 
observed for 4 min between every two or three target scans to mon- 
itor any time-dependent phase and amplitude fluctuations of the 
telescope. Flux density calibration was performed using observa- 
tions of 3C48 or 3C286, at the beginning and end of each observing 
session. The task SETJY was used to calculate 610-MHz flux den- 
sities of 29.4 and 21.1 Jy respectively, using the Astronomical Im- 
age Processing Software (AlPS) implementation of the lBaars et^ 
(1I977.) flux density scale. 

Data reduction w as performed in a simila r way to our pre- 
vious GMRT surveys (iGarn et al.|[200l l2008al lbl). Initial editing 
of the data took place separately on each sideband with stan- 
dard MPS tasks, to remove bad baselines, antennas, and channels 
which were suffering from large amounts of narrow band inter- 
ference, along with the first and last integration periods of each 
scan. The flux calibrators were used to create a bandpass correc- 
tion for each antenna. In order to create a pseudo-continuum chan- 
nel, five central frequency channels were combined together, and an 
antenna-based phase and amplitude calibration created using obser- 
vations of Jl 6 13+342. This calibration was applied to the original 
data, which was then compressed into 1 1 channels of bandwidth 
1.25 MHz, small enough that bandwidth smearing was not a prob- 
lem for our images. The first and last few spectral channels, which 
tended to be the noisiest, were omitted from the data, leading to 
an total effective bandwidth of 13.75 MHz in each sideband. Fur- 
ther flagging was performed on the 1 1 -channel data set, and the two 
sidebands combined using UVFLP jGarn et al.l2007h to improve the 
uv coverage. Baselines shorter than 1 kA were omitted from the 
imaging, since the GMRT has a large number of small baselines 
which would otherwise dominate the beam shape. 

Each pointing was divided into 31 smaller facets, arranged in 
a hexagonal grid, which were imaged individually with a separate 
assumed phase centre. The large total area imaged (a diameter of 
l.°8, compared with the full width at half-maximum of the GMRT, 
which is ~0.°74) allows bright sources well away from the observed 
phase centre to be cleaned from the images, while the faceting pro- 
cedure avoids phase errors being introduced due to the non-planar 
nature of the sky. All images were made with an elliptical restor- 
ing beam of size 6.5 x 5 arcsec^, PA +70°, with a pixel size of 
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Figure 1. The rms noise map for the GMRT ELAIS-N2 survey. The grey- 
scale ranges between 80 and 400 /iJy beam~\ and the contour levels are 
drawn at 120 and 240 /iJy beam~^ (black and white lines respectively). 



1.5 arcsec to ensure that the beam was oversampled. The images 
went through three iterations of phase self-calibration at 10, 3 and 
1 min intervals, and then a final round of self-calibration correct- 
ing both phase and amplitude errors, at 10 min intervals, with the 
overall amplitude gain held constant in order not to alter the flux 
density of sources. 

The offset in the GMRT primary beam seen in our previous 
surveys was again detected and corrected, using the method de- 
scribed in lGarn et The 13 pointings were mosaicked to- 
gether, taking into account the offset primary beam and weighting 
the final mosaic appropriately by the relative noise of each point- 
ing. The mosaic was cut off at the point where the primary beam 
correction dropped to 20 per cent of its central value, a radius of 
0.°53 from the centre of the outer pointings. 

The rms noise for the individual pointings, before correction 
for the GMRT primary beam, was between 83 and 90 /iJy beam~^ , 
compared with the expected noise level of 80 jUJy beam~^ (calcu- 
lated using Equation 1 from iGarn et 'Zl l2007h . and the rms noise 
map for the mosaic is shown in Fig.[T] 



2.2 Observations of ELAIS-Nl and the Lockman Hole 



610-M Hz GMRT observations of the ELAIS-Nl field (iGarn et all 
l2008ah and the central region of the Lockman Hole field 
jGarn et aDl2008bh have previously been presented. A brief sum- 
mary of the observations is given here; for further details, the rele- 
vant papers should be consulted. 

The ELAIS-Nl survey covers ~9 deg^ (the majority of the 
Spitzer observation region) with a resolution of 6 x 5 arcsec^, 
PA +45°. The rms noise level of the majority of the survey was 
~70 ^Jy beam~^ before primary beam correction, with a smaller 
'deep' region having ~40 /iJy beam~^ rms. 

The Lockman Hole survey covers ~5 deg^, within the central 
area of the region observed by Spitzer. The survey resolution was 
6x5 arcsec^, PA +45°, and the noise level ^60 jUJy beam~^ over 
the majority of the region, with a small area having significantly 
greater noise due to the presence of a nearby bright radio source. 



3 DATA SELECTION 

3.1 Creating an infrared sample 

The infrared dat a comes from t he band-merged SWIRE source cat- 
alogues of Sura ce et aP (l2005l), a r id the photometric redshift cat- 
alogue of Rowan-Robinson et al. (2008). In order to be included 
in the band-merged catalogue, sources were required to be de- 
tected above a signal-to-noise ratio (SNR) of 10 and 5 at 3.6 ^m 
and 4.5 fim (approximately equivalent to 10 /iJy in both bands). 
There are fewer sources present in the photometric redshift cata- 
logue than the band-merged catalogue, due to the need for sufficient 
spectral information that a redshift could be determined. Sources 
were required to be present in both the band-merged and photo- 
metric redshift catalogues, due to the need for a redshift, and the 
increased amount of information present in the band-merged cata- 
logue. This provided information on infrared flux densities in the 
four IRAC bands, along with 24-^m and 70-;Um flux densities, 
and five-band optical photometry (Ug'r'i'Z') from the Wid e Field 
Survey fWFSi'McMa hon et al.l[200ll : llrwin & Lewisll200l ') for the 
ELAIS-Nl and ELAIS-N2 fields, with four-band (Ug'r'i^) photom- 
etry for the Lockman Hole from the SWIRE photometry program. 

Sources must also be classified as 'star-forming', and have an 
estimate of the SFR. This classificat ion and estimate comes from 
the template fitting performed by Rowan-R obinson et al ] (l2008h 
during their photometric redshift calculations, and rejects some 
sources that are identified as having their infrared energetics domi- 
nated by AGN activity. An infrared detection at both 24 and 70 /im 
is also required in order to constrain the thermal emission from 
dust, and to have sufficient information on the region of the in- 
frared SED where the bulk of the bolometric energy output of star- 
forming galaxies is seen. Selecting galaxies at both 24 and 70 ^m 
increases the accurac y of the infrared bolometric lu minosity calcu- 
lation performed by Rowan-Robinson et al which would 
otherwise have a factor of ~2 uncertainty, and improves the accu- 
racy of their SFR estimates. Finally, we require the sources to lie 
within the area that was observed in one of the GMRT radio sur- 
veys. The numbers of sources in each field that satisfy the selection 
criteria are given in Table [T] 



3.2 Matching to the radio data 

In creating the public catalogues for the ELAIS-Nl and Lockman 
Hole GMRT surveys, a two-stage selection criterion with an in- 
creased SNR require ment near to brigh t sources was implemented 
(for more details see Garn et 'Zl l2008ah - this was due to problems 
encountered from residual phase errors near to the brightest radio 
sources in the survey. For this study, we have known source po- 
sitions from the infrared and optical detections, and so relax this 
criterion to a SNR of 4. Radio sources were matched to the infrared 
data if the source centres were within a distance r of each other 
- if there were potential multiple matches, the closest match was 
selected. There will be some incorrect associations between the in- 
frared and radio catalogues - we estimated this value by shifting 
the coordinates of all radio sources by 10 arcmin in declination and 
seeing how many (incorrect) associations were then made at the 
same matching radius r. We required the false matching rate to be 
less than 5 per cent, which corresponded to r = 2.5 arcsec. The num- 
ber of galaxies with a radio detection within 2.5 arcsec are listed in 
Table [TJ with a total of 510 galaxies in all three survey fields that 
satisfied the infrared and radio selection criteria. 
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Table 1. The number of sources in each sub-sample - see Sections [3ni and [T2| for further details on sample selection. 





ELAIS-Nl 


ELAIS-N2 


Lockman Hole 


Total 


Infrared band-merged catalogue (Surace et al. 2005) 


282711 


126056 


323044 


731811 


Infrared photometric redshift catalogue (Rowan-Robinson et al. 2008) 


218117 


125364 


229238 


572719 


Band-merged and photometric catalogues 


121422 


61610 


132736 


315768 


Star-forming sources only, with SFR estimate 


49160 


24831 


42517 


116508 


Star-forming, and 24-jUm detection 


11371 


5725 


12458 


29554 


Star-forming, 24-jUm and 70-/im detections 


2470 


1344 


2939 


6753 


Star-forming, 24-/im and 70-/im detections, within radio coverage area 


2204 


1317 


1891 


5412 


Star-forming, 24-/im and 70-/im detections, plus 610-MHz radio detection 


279 


63 


168 


510 



3.3 Sample selection biases 

We have rejected sources from the infrared- sele cted samples that 
were classified bv lRowan-Robinson et alJ (l2008h as having optical 
and infrared photometry that is best fitted by an AGN template; 
however there are several examples of gala xies that undergo con- 
current A GN and star formation activity fe.g. lRichards et al. l2007l : 
iNorris. M iddelberg & Boyle 2007) which will not necessarily be 
identified during this process, and this selection will not remove 
sources which have a significant contribution to their radio emis- 
sion from AGN activity, but are not clearly AGN- like at shorter 
wavelengths. We return to the question of AGN contaminants in 
Sectionim 

Selecting galaxies at both 24 and 70 jim places a signifi- 
cant restriction on the sample size, as seen in Table [T] Approxi- 
mately 25 per cent of the star-forming galaxies with an estimate 
of SFR were detected in the 24-/im band, and of these galaxies, 
approximately 25 per cent were detected at 70 jim. The infrared 
emission at 24 jim is dominated by transie ntly heated polycyclic 
arom atic hydrocarbon (PAH) molecules (e.g. Leger & PugeJl 19841 : 
IPug et & Leger 1989; Draine 2003, and references therein) and is 
strongly related to the amount of ongoing star formation, while the 
emission at 70 jim is dominated by thermal emission from large 
dust grains. Both detections are considered necessary in order to be 
sure of the accuracy of the template fitting and SFR estimation. 

Fig. 12^ shows the redshift distribution of lA-jim selected 
galaxies. The bulk of galaxies are at z < 0.5, although there is a 
secondary peak at z ~ 1 due t o the shifting of the 11.7-jUm PAH 
feature into the 24-jUm band (iRowan-Robinson et al ] l2008h . The 
70-/im flux limit for the SWIRE surveys is significantly higher than 
the 2 4-jUm limit (30 mJy compared with 0.45 mJy; Surace et al.. 
l2005h . and by selecting galaxies with 70-/im fluxes, the sample 
is biased towards galaxies that are brighter at mid-infrared wave- 
lengths, and have high SFR. The secondary redshift peak at z ~ 1 
is reduced in strength for the sample selected at lA-jim and 70- 
/im, but remains prominent - this is due to the global decrease in 
SFRD of the Universe, seen between z ^ 2 and the current day (e . g . 
iLiUv et al.l[T993 : iMadau etaD llQQ^: 'Hopkins & BeacomI 12006'). a 
consequence of which is that the highest SFR galaxies are likely to 
be found at higher redshift. The SFR distribution of the two samples 
is given in Fig.[2l3, and shows that a greater fraction of galaxies in 
the 24-^m and 70-;Um selected sample have SFR > 100 Mq yr~^ 
than when only a 24-/im detection is required. Fig. [3] illustrates 
that high-SFR sources are more likely to be found at higher z in our 
sample. The upper left of the diagram is relatively unpopulated due 
to the lower space density of high SFR sources in the local Uni- 
verse compared to at z ~ 2, while the lower right of the diagram is 
unpopulated principally because of the flux limit of the 70-^m data 
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Figure 2. The fractional distribution of (a) redshift and (b) SFR for galaxies 
selected at 24 /im (soHd line), 24 and 70 /im (dashed line), and 24 jUm, 
70 /im and 610 MHz (dotted line) - see Section [331 for further details. 

used in the iRowan-Robinson et ID (l2008h catalogue. The curved 
line indicates the minimum detectable SFR at different redshifts, 
based upon the 70- /im flux only, using the relationship between a 
monochromatic MIPS flux density and SFR given by Rieke et al.l 
(l2009h . There are a few sources seen below this limit, due to the 
method used by iRowan-Robinson et al.l (l2008h to estimate SFRs 
(see Section 14.3.31) , which considers more information than just a 
single 70-/im detection. Note that the shape of the SFR limit de- 
pends both on the intrinsic limiting luminosity for a flux-limited 
sample, and the /^-correction caused by the slope of the mid-infrared 
SED for a typical starburst galaxy. 

The photometric redshift catalogue o fE owan-Robinson et al. 
(2008) contains a considerable number of sources with an inferred 
SFR of > 1000 Mq yr~^. Such extreme SFRs are unusual, and 
for man y sources are likely to be unphysical, but are not unknown 
(see, e.g. ICapak et al.l2008l . for a discussion of a spectroscopically- 
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Figure 3. The relationship between z and SFR for the 510 galaxies detected 
at 24 jim, 70 jim and 610 MHz. The solid line represents the limiting SFR 
detectable at different redshifts, using the relationship between monochro- 
matic 70-jUm flux density and SFR given bv .Rieke et al.. 12009.) . 



confirmed source with SFR - 4000 M© yr-^). We note that SFRs 
derived from either 24-/xm or lO-jim fluxes are consistent with 
these high values for most of the extreme SFR sources in our sam- 
ple (se e Section |4. 3. 3 K although the greatest SFR that lRieke et alJ 
(l2009h model is - 1600 M© yr-^ and inferred SFRs above this 
value will be very poorly constrained. We also note that SFRs de- 
rived independently from radio data are consistent with these high 
values (see Section |431 ), although in both cases an error in the cal- 
culated photometric redshift may be leading to a significant over- 
estimation of the SFR. The small number of sources in our sample 
with a very high inferred SFR do not affect our later conclusions 
on variations in the relationship between radio luminosity and SFR 
with redshift. 

Selecting galaxies with photometric redshifts, rather than re- 
quiring spectroscopic redshifts, improves the homogeneity of the 
sample - all galaxies with sufficient photometry will be included 
in the study, rather than limiting the sample to galaxies that are 
sufficiently bright at one wavelength that they were selected for 
spectroscopic followup. Photometric redshifts are inherently less 
accurate than spectroscopic equivalents - however, there were con- 
siderable numbers of spectroscopic redshifts available for valida- 
tion in the SWIRE fields, and Rowan-Robinson et al. (2008) esti- 
mate an typical error of 4 per cent in (zphot - ^spec)/(l +^spec), 
with ~1 per cent of the photometric redshifts being catastroph- 
ically incorrect. Where both spectroscopic and photometric red- 
shifts are available, we use the spectroscopically determined value, 
and we return to the possibility of catastrophic redshift errors in 
Section l433l 

The radio surveys used in this work are much less sensitive to 
star-forming galaxies than the infrared surveys, as can be seen by 
the low number of source detections in the bottom row of Table [T] 
The percentage of infrared- selected galaxies that are detected at 
610 MHz is 4.8 per cent (ELAIS-N2), 8.9 per cent (Lockman Hole) 
and 12.7 per cent (ELAIS-Nl), with the fields that have lower 610- 
MHz noise levels having a greater percentage of source matches. 
Due to the po sitive correlation that is seen bet ween radio luminosity 
and SFR (e.g. lCondon & Yinll99d : lBelll2003u and see Sectionl42l), 
additionally requiring a detection at 610 MHz will bias the sample 
against galaxies with very low SFR. This can be seen in Fig. [Jt - 
the fraction of galaxies with SFR < 1 Mq yr~^ is significantly re- 
duced compared with the infrared-only samples. A study of sources 
below the radio detection limit, r equiring a statistical process such 
as me dian source stacking (e.g. iBovle et al.l [20071 : iBeswick et al.l 
l2008h . is beyond the scope of this paper, but has been discussed 
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Figure 4. Variation in (a) and (b) q'-jQ with redshift, along with the local 
median value of (black line with error bars). Sources which have been 
identified as being AGN-like from their deviation below the median value of 
^24 are indicated by crosses, '+' - see Section lTTI for further details. Open 
circles, 'o', mark the sources that have been identified as being radio-quiet 
compared with their SFR - see Section l43] for further details. 



in a separate work (iGarn & Alexandei]l2009h . We will return to the 
galaxies where a non-detection in the radio is significant in Sec- 
tionl431 



4 RESULTS 

4.1 Identifying potential AGN contaminants 

There is a tight relationship between the infrared and radio 
luminosities of galaxies that are predominantly powered by 
star-formation processes (the 'infrared / radio coirelatio n' ; e.g. 
iHelou et al.|[T985l : ICondon et aLlll99ll : iMurphv et al.ll2006h . which 
is not f ollowed by sources that have significant a mounts of AGN 
activity (Sopp & Alexander "l 99 iIi IRov et al ."l998V The correlation 
has previous ly been used as a discriminant between the two types of 
source (e.g. Ilbar et al.ll2008l : [Sevmour et aLll2008h . and in order to 
use this technique on our sample, we consider the monochromatic 
parameter given by 



^{r = log 



10 



>^1R 
*^610 



(1) 



where 5ir is the infrared flux density in either the 24-/im or 70- 
^m bands, and S^iq is the radio flux density at 610 MHz. Note 
that this definition is slightly different from previous works (e.g. 
lApDleton e t"aDl2004h which typically use 1.4-GHz flux densities. 
Fig.|4]shows the variation in and q^Q with redshift for the galax- 
ies in this sample. The dispersion in the individual values is 
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greater than typically seen at 1.4 GHz (lAppleton et al.ll2004h . The 
principal reason for this is that for galaxies which have both AGN 
and star-forming activity, the 610-MHz flux density mainly stems 
from the nuclear source, while the 1. 4-GHz flux density come s 
mainly from star formation activity dMagliocchetti et all l2008l) . 
AGNs were identified by find ing sources that were a better fit to 
the selected AGN templates of iRowan-Robinson et alJ (l2008h than 
to elliptical galaxy templates, as discussed in Section[3j these were 
then removed from the sample. Any sources that have their radio lu- 
minosity dominated by an AGN, but do not clearly show up as hav- 
ing infrared AGN templates, will not have been removed through 
this method. These 'missed' AGN sources will be radio bright, and 
therefore have a low value of q'^^. 

We have not converted the data in Fig.|4lto a consistent rest- 
frame wavelength (a '^-correction'), and so the observed value of 
q'^^ will vary with the redshift of the galaxy, as different portions 
of the redshifted SED enter into our observation bands. In order 
to quantify how the typical value of q'^^ varies, we calculate the 
median value of q'^^ in ten logarithmically- spaced redshift bins be- 
tween z = and 2, and show this as the solid line in Fig. IH along 
with error bars indicating the local standard deviation o. The me- 
dian value of ^jj^ shows a steady decline with redshift at both 24 
and 70 jim, due to the shape of a typical galaxy SED, which has 
a steeper spectrum in the mid-infrared than in the radio. We define 
sources which have a value of below the local median value by 
more than 2(7 to be 'AGN-like' - the typical value of o for a bin 
is 0.4, making these sources a factor of ~ 6 mor e radio-bright tha n 
is typical. This criterion is similar to that used bv llbar et"ar (l2008h . 
although they used lo below the mean value of ^24 as their cut-off 
value. The median estimator is more robust to the presence of out- 
liers than the mean; since it is these outliers which we are trying 
to identify, the median is a more appropriate statistic to use. There 
are 25 sources out of 510 (5 per cent) which satisfy this criterion, 
seven of which are more than 3(7 away from the local median. The 
AGN-like sources are identified in Fig.|4l For comparison, there are 
only five sources that lie above the median value of by more 
than 2(7, and none by more than 3(7. A similar selection for AGN- 
like sources could be made in terms of differences from the median 
value of - however, most of the sources that show significant 
deviations at 70 /im have already been identified through their ^24 
deviation, and since the resolution of the lA-jim data is a much bet- 
ter match to the 610-MHz data we choose to consider deviations 
from ^24 only. 

All of the AGN-like sources were visually inspected in the ra- 
dio images, in order to look for extended jet/lobe structure. One 
source was potentially the central nucleus of an extended radio 
source, but none of the other sources showed evidence for any 
extended radio structure. Fig. [5] shows optical, IRAC and MIPS 
50 X 50 arcsec^ images centred on each of the AGN-like sources. 
The sources are typically compact, although a few (e.g. the fourth 
source on the top row of Fig. [5]) show signs of extended emission 
in the optical images. Throughout the rest of this work, the AGN- 
like sources will be retained in the sample, but identified as poten- 
tial contaminants in order to compare them to the remainder of the 
population. 

4.2 Star formation rate estimates 

The SFR estimates for the galaxies in our sample come 
from the SED fitting and photometric redshift estimation of 
iRowan-Robinson et al.l (l2008h . A further method of estimating the 
SFR, ^, from the 1.4-GHz luminosity of galaxies is described by 
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Figure 6. The relationship between ^-corrected 610-M Hz radio lum inosity 
and SFR (grey diagonal crosses) and the prediction of 'Bell' (EoO?) shifted 
to 610 MHz assuming a = 0.8 (solid line). The sources identified as AGN- 
like are shown as crosses, '-1-', and the sources identified as being radio-quiet 
(see Section |43^ are shown as open circles, 'o'. Significant non-detections 
are indicated by upper limits - see Section |43] for further details. 



iBei] (l2003h . calibrated from the total infrared SFR for galaxies 
with L > L* (defined as having an infrared luminosity Ljr ~ 2 x 
lO^L©). Assuming a radio spectral index a = 0.8 (e.g. iCondonI 
[1992), where we define a such that the flux density S = SqV we 
convert this relationship to a 610-MHz equivalent: 



: 2.84x10-22 (-^ 



WHz- 



(2) 



Moyr-i 

for Lgio > Lc (where Lc = 3.3 x 10^^ W Hz~^ is the luminosity at 
610 MHz of a ~ L* galaxy, with ^ ^ 1 Mq yr"^), and 

^ \ 2.84x10-22 f L610 



Moyr- 



0.1+0.9(L6io/Lcr3 VWHz 



(3) 



forL6io <^c. 

Fig. [6] shows the relationship between 610-MHz radio lu- 
minosity and SFR. The radio data were ^-corrected from an 
observation-frame flux density into a rest- frame luminosity using 



(4) 



where we assume that a = 0.8, and where is the luminosity- 
distance corresponding to redshift z. The relationship given by 
Equations |2] and [3] is overlaid on the data, and shows a very good 
agreement with the overall trend for sources which are detected 
in the radio images. Because the sample has few galaxies with 
SFR less than 1 Mq yr"^, we are unable to confirm the existence 
of a break dX L < Lc. The sources that were previously identified 
as AGN-like all have a higher radio luminosity than predicted by 
Equation |2] (which applies to star-forming galaxies only), although 
some do not deviate far from the overall trend. 

We define the ratio of Lgio/^ to be the 'specific ra- 
dio luminosity' of a galaxy, and calculate the median value of 
logio(L6io/^), for all sources with radio detections, and for 
sources with SFR > 1 Mq yx~^ which are not classified as being 
AGN-like. The medians and standard deviations are 21.51 ±0.61 
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Figure 5. Optical and infrared images of the 25 AGN-like sources, retrieved from the NASA/IPAC infrared science archive. All images are 50 x 50 arcsec^ in 
size. Top row: /-band images. Middle row: 3.6-/im images. Bottom row: 24-/im images. 
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Figure 8. The minimum SFR necessary for sources to be identified as radio- 
quiet, based upon 610-MHz noise levels of 40 /iJy beam~^ (solid line) and 
90 /iJy beam~^ (dashed line). The vertical dotted line indicates a redshift 
of 0.2. 



and 21.48 ±0.55 respectively, consistent with the iBelll (l2003h 
value (shifted to 610 MHz) from Equation [2l of 21. 55. In contrast, 
the ra dio-luminosity / SFR relationship given bv ICondon &^ 
(Il990h . which is based around the ratio of the non-thermal Galactic 
luminosity at 408 MHz and an estimate of the Galactic supernova 
rate, predicts a value of 21.17 at 610 MHz - also within one stan- 
dard deviation of our results, but a less good fit to the observations. 



4.3 Radio- quiet sources 

There are several sources which are detected in the radio images, 
but show significantly lower values of radio luminosity in Fig. [6l 
than would be expected from their SFR and Equation [2l We de- 
fine 'radio-quiet' sources as having a value of log lo (1^6 10 /^) that 
is at least 2(7 below the median value, calculated from non-AGN 
sources only. There are 20 sources satisfying this criterion, which 
are marked as such on Fig. [6l The sources are at low redshift (see 
Fig. 131), with 90 per cent being at < z < 0.2. 

As discussed in Section [331 the Spitzer surveys are more sen- 



sitive to star-forming galaxies than the radio surveys - of the 5412 
infrared- selected sources, only 510 were detected at 610 MHz. 
In order to identify which of the non-detections are significant, 
we determine upper limits on 610-MHz flux density for each 
source, based upon four times the local noise level, and place up- 
per limits on the radio luminosity using Equation |4l We then de- 
termine an upper limit on the value of logio(L6io/^) using the 
iRowan-Robinson et al.l (12008') SFR values - for the majority of 
sources, this value lies above or within the distribution of detected 
sources shown in Fig. [6l (i.e. the radio surveys are not deep enough 
to detect the source if it has a 'typical' specific radio luminosity); 
we consider these non-detections to be insignificant, and do not dis- 
cuss them further. There remain a small number of sources which 
have an upper limit on specific radio luminosity that is more than 
2(7 below the value predicted by Equation |2l (i.e. we would expect 
to detect these sources if they had a 'typical' specific radio lumi- 
nosity). There are 48 significant non-detections, which are shown 
on Fig.[6l and appear to be similar to the 'radio-quiet' sources pre- 
viously identified. 

There are a few potential explanations for these sources, which 
we consider in turn: 

(i) we have significantly under-estimated the radio luminosity 
for these sources fSection B.3.11 ): 

(ii) the sources are genuinely radio-quiet, either due to a recent 
burst of star formation, or due to a suppression of the radio lumi- 
nosity within these galaxies (Section 14.3.21 ): 

(iii) the SFR has been significantly over-estimated by 
iRowan-Robinson et all (l2008h (Section 14331 ). 



4.3.1 Under- estimation of the radio luminosity 

Any significant under-estimation of the radio flux density would 
lead to a lower calculated value of the specific radio luminosity of 
galaxies. We inspected all of the radio-quiet sources in the radio 
images in order to test whether the radio counterparts were close 
to bright sources, or in noisy parts of the image, but no system- 
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Figure 7. Optical and infrared images of the 20 radio-quiet sources, retrieved from the NASA/IPAC infrared science archive. All images are 50 x 50 arcsec^ 
in size. Top row: /-band images. Middle row: 3.6-/im images. Bottom row: 24-/im images. 



atic effects were apparent. The sources are mainly at low redshift 
(90 per cent of the radio-quiet galaxies, and 94 per cent of the sig- 
nificant non-detections are at < z < 0.2), and galaxies are likely 
to have a greater angular size in the local Universe than at higher 
redshift. Fig. [7] shows that some of the radio-quiet sources appear 
to be extended in the optical and infrared images, suggesting that 
there may be diffuse radio emission near to or below the noise level 
which has not been included into the measured value of flux den- 
sity. 

The low redshift of these sources can be explained purely 
through selection effects - the SFR that is required for a source to 
be identified as either radio-quiet or as a significant non-detection 
depends on both the redshift of the source and the local noise level 
of the radio image, as illustrated by Fig. [8] At z = 0.2, a galaxy 
needs to have a SFR of at least 60 Mq yr~^ to be classified as 
radio-quiet in the most sensitive portion of the ELAIS-Nl image, 
with noise level of 40 jUJy beam~^, or alternately a SFR greater 
than 135 Mq yr~^ for a noise level of 90 ^Jy beam~^, more char- 
acteristic of the ELAIS-N2 survey. Fig. [Jt shows that most of the 
24-^m and 70-^m selected sources have SFRs below these values, 
and so could not be identified as being significantly radio-quiet. 

Even if the host galaxy is extended, the majority of star for- 
mati on in h igh SFR galaxies tends to come from a nuclear starburst 
(e.g. lKennicuttlll99 8), and the bulk of the radio emission would not 
be expected to be spread out over a large area. In order to test for 
extended radio emission we measured the flux density of each of 
these sources within a circular aperture centred on their known lo- 
cation, using the technique described in Garn & Alexander ( 2009), 
and calculated the ratio of flux density measured within the aper- 
ture to the flux density listed in the catalogue. For the full sample of 
sources, the flux density measured within apertures of varying sizes 
was comparable to the catalogued flux density, from which we con- 
clude that most of the sources have little extended radio emission. 
The AGN-like sources also follow this trend, but the radio-quiet 
galaxies show evidence for a slight increase in flux density within 
apertures of > 15 arcsec radius, implying the existence of some dif- 
fuse radio emission which has not been catalogued. The increase is 



less than a factor of two; this is not a large enough effect to make 
these galaxies significantly radio-quiet. 

The conversion between flux density and luminosity relies 
upon an assumed radio spectral index of 0.8 for all sources; without 
multi-frequency radio data it is not possible to estimate this value 
more accuratelv. lOarn et an([2007 ) demonstrate that the spectral in- 
dex distribution of galaxies detected at both 610 MHz and 1.4 GHz 
peaks around a = 0.8, and various studies fe.g. lBondi et al.ll2007l : 
llbar et allliooi) have concluded that there is no significant varia- 
tion in the mean spectral index of radio sources at flux densities 
down to ~ 100 /iJy. Even for an extreme error in the spectral index 
of (7a = 0.5, for a galaxy at z = 2, the inferred luminosity would 
not change by more than a factor of two; again, this is not a large 
enough effect to explain the radio-quietness of these sources. 

4.3.2 Physical effects 

There are a few potential reasons why galaxies could be radio-quiet 
with respect to their SFR; recent star-burst activity, or suppression 
of the radio luminosity due to a low confinement of cosmic rays, 
potentially caused by low magnetic field strengths. If a galaxy had 
just begun a burst of star formation activity that has enhanced the 
SFR, but not yet enhanced the supernova rate, it would appear to 
be radio-quiet - since the supernova rate lags the SFR by ~ 30 Myr 
(the approximate lifetime of 8 Mq stars), this would mean that 
sources were being observed very close to the beginning of the star- 
burst period. Some of the sources (e.g. the second source in row two 
of Fig.|7]) show signs of potential on-going mergers in their optical 
images, which would be an indication of triggered star formation, 
although this is not the case for all sources. A full analysis of the 
possibility of mergers in the sample is beyond the scope of this pa- 
per, although it seems unlikely that all of the radio-quiet sources in 
our sample could be undergoing merger activity. 

The best- studied star burst galaxy in the local U niverse is 
Arp 220, at z = 0.018126 dde Vaucouleurs et aLl[T99lh with SFR 
of 240 Mq yr~^ , calculated from it s FIR luminosity dSanders et al.l 
I2003D and the relationship given in iKennicutd (Il998h . It has a 1.4- 
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GHz flux density of 0.326 Jv dCondon et alJl2002h and a 365-MHz 
flux density of 0.435 mJy jPoudas et al.lll996h . giving it a spectral 
index a ^ 0.2. This equates to a luminosity of 2.5 x 10^^ W Hz~^ 
at 610 MHz, and a ratio of log^g (^6io/^) of ^ 21, within la 
of the median value found in Section 14.21 thus Arp 220 would 
not have been identified as being radio-quiet through this diagnos- 
tic. However, this should not be surprising - A rp 220 has a large 
number o f compact radio sources visible fe.g. Smith et all 1 19981 : 
iLonsdale et al. 2006) which are thought to be radio supernovae 
formed due to the starburst activity - this system has therefore al- 
ready passed th e phase at whi ch it could be identified as radio-quiet. 
The findings of lWilson"et aD(2006) agree with this interpretation - 
there are a number of star clusters with ages up to around 500 Myr 
within Arp 220, suggesting that starburst activity started a long time 
ago. Any starburst galaxies identified through through this diagnos- 
tic would have to be significantly earlier on in their starburst phase 
than Ar p 220. 

Carilli et al ] (I2OO8I1 have identified evidence for a possible 
decrease in the conversion factor between radio luminosity and 
SFR for high redshift (z ~ 3) Lyman break galaxies, based upon 
a stacking analysis of sources in the COSMOS field (although 
see Garn & Alexander 2009, for a discussion of the biases present 
in radio stacking experiments). Their favoured explanation of in- 
creased relativistic electron cooling due to inverse Compton scatter- 
ing off the cosmic microwave background (CMB) would not apply 
to our low-redshift radio-quiet galaxies, since the energy density 
in the CMB increases with redshift as (1 +z)^ and is compara- 
tively unimportant in the local Universe. A potential explanation 
for low-redshift radio-quiet sources would be that they have lower 
magnetic field strengths than is typical, leading to a decreased 
electron confinement (for the optically-thin scenario discussed by 
IchT & Wolfendale 1990) and a lower fraction of the total energy in 
cosmic rays being radiated within a galaxy compared with a source 
having a greater 5-field. However, the radio-quiet sources follow 
the same infrared-radio correlation as the main sample (shown by 
Fig. |4]) and a decrease in radio luminosity would therefore also re- 
quire a corresponding decrease in the infrared luminosity; if this 
were the case, the measured SFR (which was calculated from the 
infrared luminosity) should also decrease, and sources would not 
be identified as being radio-quiet. 



4.3.3 Over-estimation of the star formation rate 

The SFR for these galaxies has been estimated u sing a complex 
method, described in lRowan-Robinson et al.l(l2008h : 

(i) optical and IRAC fluxes are used to find the best-fitting tem- 
plate spectrum for each galaxy, from a selection of six templates, 
and photometric redshifts are calculated; 

(ii) sources with an 'infrared excess' at A > 8 fim have their 
bolometric infrared luminosity Ljr estimated from the long wave- 
length data; 

(iii) 60-/im luminosities, Lgo, are calculated from Ljr, using 
bolometric correction factors of 3.48, 1.67 and 1.43 for cirrus, M82 
and Arp 220 galaxy templates respectively. 

(iv) SFRs are then estimated from the 60-jUm luminosity, using 
the conversion factor of 

where £ describes the fraction of UV light absorbed by dust, and is 
taken to be 2/3. 




Figure 9. Comparison of spectroscopic and photometric redshifts for all 
galaxies in our 24-/im and 70-/im selected sample with spectroscopy (grey 
diagonal crosses), the radio-quiet sources (open circles) and the significant 
non-detections (stars). The 1:1 line is indicated for reference. 



iRowan-Robinson et "Zl (l2008h estimate an uncertainty in the 
total infrared luminosity at z = 0.2 (the approximate redshift of the 
radio-quiet sources) of 0.2 dex - this is too small to be responsi- 
ble for the deviations which are seen. A more plausible scenario is 
for source confusion to play an important role - the resolution at 
70 jUm is only ~ 18 arcsec, much poorer than in the optical, IRAC 
or 24-/im bands. If there are two distinct sources present in the op- 
tical images, only one of which has been assigned a far-infrared 
counterpart (with the combined infrared flux from both) then that 
source will have its SFR overestimated in the photometric redshift 
catalogue, while the other source will not appear in the sample due 
to its lack of an infrared detection. Inspection of Fig. |7] shows sev- 
eral sources where confusion could potentially be having an effect 
on the SFR estimation. However, it could only plausibly increase 
the flux density by a factor of ~ 2, and is not likely to lead to an 
order of magnitude discrepancy. 

The conversion between Ljr and L^q has uncertainties in it, but 
again there is only a factor of ~ 2 difference between the three tem- 
plates that are used. Likewise, the parameter e used to calculate the 
amount of young star formation coul d be in error for galaxies with 
a significant old stellar population (Bell""2003|; iRowan-RobinsonI 
2003) which would lead to an over-estimation of the SFR, but nei- 
ther of these factors could explain such a large difference in SFR 
for these sources. 

The fact that the radio-quiet sources follow the same infrared 
/ radio correlation as the main sample, but yet show a disagree- 
ment between their radio luminosity and in frared-derived SFR im- 

i )lies t hat the template-fitting carried out bv lRowan-Robinson et al.l 
20081) to estimate photometric redshifts and total infrared luminos- 
ity may have been unsuccessful for a small number of sources - if 
this were the case, then the values of SFR derived from this process 
would be inaccurate. In order to test for this, in Fig.|9]we compare 
the photometric and spectroscopic redshifts for all of the galax- 
ies which have available spectroscopy, and were selected at 24 and 
70 jUm. All 20 of the radio-quiet sources, and 46 of the 48 signifi- 
cant non-detections have spectroscopic redshifts, in contrast to the 
fraction of spectroscopic red shifts available for the whole s ample 
(576/5412; 11 per cent). The iRowan-Robinson et aP (l2008h crite- 
rion for defining a catastrophic redshift error is Alogio(l +z) = 
±0.06, or 15 per cent; of the 66 radio-quiet or significant non- 
detected sources with spectroscopic redshifts, 30 (45 per cent) were 
outside this threshold. In contrast, only four of the sources which 
were detected at 610-MHz and not found to be radio-quiet were 
classified as having a catastrophic redshift error. We exclude these 
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Fi gure 10. The relationship betw een SFRs estimated by the template-fitting 
of Rowan-Robinson et al H2008t), and SFRs estimated from the relation- 
ships eiven in lRiek^^] ( 2009), using (a) 24-;Um flux densities; (b) 70- 
}im flux densities. Galaxies in the 24-/im and 70-/im-selected sample are 
shown as grey dots, with the radio-quiet sources indicated by open circles, 
and the significant non-detections shown as stars. 



34 sources from the remainder of the analysis, as the SFRs calcu- 
lated using inaccurate templates will be unreliable. 

In order to test whether the SFRs for the remainder of the 
radio-quiet and significant non-detected sources are reliable, we use 
relationships provided by Rieke et al. ( 2009) to calculate SFRs di- 
rectly from monochromatic flux densities measured by MIPS at ei- 
ther 24 or 70 ^m. While SFRs based around a single flux measure- 
ment are likely to be less accurate than SFRs based around bolo- 
metric luminosities, Rieke et al. (2009) believe that they should 
be accurate to within about 0.2 dex, and these estimates can 
be used as an independent verification of the accuracy of the 
iRowan-Robinson et alj (l2008 h SFRs. In Fig.[Tolwe compare these 
three SFR estimates for our lA-jim and 70-/im selected sample. 
There is some scatter in the data, but it is clear that the galax- 
ies which are radio-quiet or significant non-detections are typically 
outliers, wit h their monochromatic SFRs be ing smaller than the es- 
timates from ^ Rowan-Robinson et alj (l2008l) . 

We favour an over-estimation of the SFR in a small fraction 
of our sample for explaining the apparent low specific radio lu- 
minosity of a few of the sources, although we can not rule out a 
contribution from the other possibilities. We note that if SFRs cal- 
culated from either monochromati c 2A-\im or 70- )Um flux densities 
were used instead of the Rowan-Robinson et 3(2008) values for 
these sources, only two of the significant non-detections would re- 
main significant, and none of the radio-quiet sources would still 
be classified as such. For completeness the radio-quiet sources and 
sources with a significant non-detection, which have not been re- 



jected for poor template-fitting, will be retained on future plots for 
comparison with the rest of the population. 



5 DISCUSSION 

iBelll (l2003h presented a relationship between the SFR and 1.4-GHz 
non-thermal radio luminosity of galaxies in the local Universe. This 
was calibrated from local galaxies, assuming a direct proportional- 
ity between non-thermal radio emission and SFR for galaxies with 
LiR > L*. By constructing a sample of galaxies within the northern 
SWIRE regions, we have demonstrated that this relationship can 
be successfully applied to galaxies that are detected at 610 MHz. 
We compared the 610-MHz radio luminosity (^-corrected assum- 
ing a synchrotron radio spectrum with a = 0.8) to the SFR in 
order to calculate a parameter, L^i^/m, which we call the 'spe- 
cific radio luminosity' of star-forming galaxies. This value s hows 
good agreement with the relationship given bv lid] (l2003h . and 
a less good niatch to the more commonly-used relationship of 
ICondon & Yh] (Il990l ). We find that approximately 4 per cent of 
the detected sources do not follow this relationship, and note the 
existence of several more galaxies that would be expected to be de- 
tected at 610 MHz, but also appear radio-quiet. Of the total number 
of lA-jim and 70-/im selected galaxies which could potentially be 
detected at 610 MHz, we find 12 per cent (68/558) are significantly 
more radio-quiet than the iBelll (12003 ) relationship would suggest 
from their SFRs. We suggest a number of potential explanations 
for these sources, and favour an inappropriate conversion between 
infrared flux densities and SFRs for a small number of galaxies, 
due to an uncertainty in the photometric redshift calculations, rather 
than a real physical effect. 

Various studies appear to indicate that there has been li ttle evo- 
lution in the infrared / radio correlation at h igh redshift (^e.g. lOarretl 
l2002l : lApDleton et al Jl2004l : llbar et alJl200 8). suggesting that there 
has been little change in the physical processes which link star for- 
mation, thermal dust emission and radio luminosity, although these 
studies work with sources selected in the radio and are therefore 
biased against radio-quiet galaxies. In this work we do not directly 
test for variations in the luminosity correlation with redshift. How- 
ever, the SFR estimates provided by Rowan-Robinson et al. ( 20081) 
allow a probe of any evolution in the physical link between SFR 
and radio synchrotron luminosity at high redshift, which has the 
advantage over studies which look at either the luminosity or flux 
density correlation in that we are working with one physical param- 
eter (the SFR) and one potentially variable process (synchrotron 
radiation), rather than two processes (thermal dust emission at a 
given wavelength and synchrotron radiation) that could both vary 
with redshift. 

In Fig.[TT]we show the specific radio luminosity of the galax- 
ies in our sample, as a function of redshift. Any variation seen in 
specific radio luminosity with z would suggest a change in one or 
more of the properties that controls the synchrotron radiation from 
star-forming galaxies - specifically, the galactic magnetic field 
strength, or the amount of confinement of the high-energy elec- 
trons. However, no such variation is seen - the majority of galaxies 
w hich a r e dete cted in the radio follow the global relationship given 
bvlidil (12003*) out to a redshift of at least 2 (a look-back time of 
10.3 Gyr) and potentially to z ~ 3. The AGN-like sources discussed 
in Section|4T|all have a high specific radio luminosity (which is to 
be expected if some of their radio luminosity is not related to star 
formation activity), but their specific radio luminosity shows no ev- 
idence for redshift evolution. The radio-quiet sources discussed in 
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Figure 11. The relationship between 610-MHz radio luminosity and SFR, 
as a function of reds hift - symbo ls are the same as for Fig. |6] The local 
relationship given by iBelll l2003h is indicated by the solid line - there is 
no evidence for variation in the radio luminosity / SFR relationship out to 
z = 2, for radio-selected star-forming galaxies. 
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Figure 12. The relationship between 610-MHz radio luminosity and 
infrared-derived SFR, as a funct ion o f SFR - symbols are the same as for 
Fig- [SI The relationship given bv lid] ( [2003b is indicated by the solid line. 



Section [431 all lie at low redshift and do not affect this conclusion, 
although we caution that due to selection effects, we are unable to 
test for the existence or the number of radio-quiet sources which 
may exist at redshifts beyond 0.2. 

As discussed in Section [TJ there is much interest in trying to 
calculate the star formation history of the Universe. Radio obser- 
vations at ~ 1 GHz are unaffected by dust, and should allow a 
reliable estimation of the variation in SFRD with z, without the 
uncertainties that are found at optical and infrared wavelengths. 
The re have been a few atte mpts to carry out this calculation (e.g . 
Mobasher et al. 1999; Haars ma etal][200Ql : [Sevmour et aP [20081 : 
Smolcic et al.ll2009 ) but all of these rely upon the untested assump- 
tion that the local relationship between radio luminosity and SFR 
can be applied successfully at higher redshift. The work presented 
here demonstrates that this assumption is valid - the radio lumi- 
nosity of the majority of radio- selected star-forming galaxies does 
track the SFR out to at least z = 2 (with some potential evidence of 
an extension out to z ~ 3 and beyond; see Fig.[TT[). 

Equally importantly, the relationship between luminosity and 
SFR has been implicitly assumed to be the same for galaxies with 
all values of SFR. Fig. [12] shows that, at least over the range 
1-10^ Mq yr~^, no significant variation in this relationship is 
seen. Galaxies undergoing extreme starbursts, with SFRs up to 
10^ Mq yr~^, and galaxies undergoing much lower star-formation 
activity with SFR of 1 Mq yr~^ show the same relationship be- 



tw een their syn chrotron emission and SFR. The deviation predicted 
bv [Bd]([2003[) below ~ 1 Mq yr ^ cannot be tested in this work, 
due to a lack of low-luminosity sources in the sample. 

We have demonstrated that the Bell (2003) relationship be- 
tween 1.4-GHz luminosity and SFR (for L > L* galaxies) can be 
used successfully on galaxies at < z < 2, although it is necessary 
to first apply a ^-correction to the observed flux density of a galaxy 
in order to convert to a rest-frame 1.4-GHz luminosity. Equation [6] 
gives a gener alised versio n of this relationship, obtained through 
combining the [Bell[ ([2003[) equation with Equation [J] 



Moyr- 



0.066 fdL{z)y f {l^z)v Y fSy 

1+z \ Mpc J vi-4GHzy VJy 



(6) 



This equation can be used to predict the SFR of a galaxy directly 
from observational radio data. 



6 CONCLUSIONS 

We have presented new 610-MHz GMRT observations of the 
ELAIS-N2 region, covering 6 deg^ with 13 pointings. The typi- 
cal noise level of each pointing is ~90 /iJy beam~^ before pri- 
mary beam correction. These radio observations have been com- 
bined with our previous surveys of the ELAIS-Nl and Lock- 
man Hole regions, along with optical and infrared surveys, in 
order to create a sample of 510 galaxies in the three northern 
SWI RE fields, which have photome tric redshifts and SFR estimates 
from [Rowan-Robinson et al.[ ([2008[) . infrared detections at 24 and 
70 /im, and radio detections at 610 MHz. 

We have used the logarithmic ratio of infrared to radio flux 
densities, ^(j^, to discriminate between sources which are powered 
by star formation processes, and sources which may have a sig- 
nificant fraction of their radio emission resulting from AGN activ- 
ity. By comparing the rest-frame 610-MHz radio luminosity to the 
Rowan-Robinson et al. ( 2008) SFR, we have identified a further set 
of sources which appear to be radio-quiet; we have demonstrated 
that these are likely to be galaxies where the template-fitting has 
been unsuccessful, leading to an over-estimation of the photometric 
redshifts and SFRs, although we discuss a number of real physical 
effects which could also be contributing to a reduction in the radio 
luminosity, such as the rece nt onset of starburst ac tivity. 

We have compared the [Condon & Yin[ ([1990*) and'EdJ ([2003[) 
relat ionships between SFR and radio luminosity, and found that the 
|Bd3 12003) prescription is a better fit to our sample. By considering 
the ratio of radio luminosity and SFR, the 'specific radio luminos- 
ity', we have tested the validity of two of the assumptions that are 
commonly made when calculating the SFR history of the Universe 
from deep radio observations, namely that a relationship calibrated 
from local galaxies is applicable to galaxies at much higher red- 
shift, and undergoing much greater starburst activity than seen in 
the local Universe. We find no redshift dependence for the spe- 
cific radio luminosity, suggesting that there has been little change 
in the physical processes linking synchrotron radiation and star for- 
mation in star-forming galaxies over the redshift range < z < 2. 
No variation in the specific radio luminosity of galaxies with their 
SFR was found for a SFR in the range of 1 to 10^ Mq yr~^. This 
implies that the link between SFR and synchrotron luminosity is 
the same for massive starburst galaxies as it is for the more q uies- 
cent g alaxies seen in the local Universe. We conclude that the [Bell[ 
([2003[) relationship between radio luminosity and SFR, calibrated 
from local galaxies, can successfully be applied to high-redshift, 
high-SFR galaxies, and present a generalised equation to link the 
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SFR of a galaxy to the observed flux density, redshift and radio 
spectral index. 
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